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The dynamical stability of W Ursae Majoris-type systems 
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ABSTRACT 

Theoretical study indicates that a contact binary system would merge into a rapidly 
rotating single star due to tidal instability when the spin angular momentum of the 
system is more than a third of its orbital angular momentum. Assuming that W UMa 
contact binary systems rigorously comply with the Roche geometry and the dynamical 
stability limit is at a contact degree of about 70%, we obtain that W UMa systems 
might suffer Darwin's instability when their mass ratios are in a region of about 0.076- 
0.078 and merge into the fast-rotating stars. This suggests that the W UMa systems 
with mass ratio q < 0.076 can not be observed. Meanwhile, we find that the observed 
W UMa systems with a mass ratio of about 0.077, corresponding to a contact degree of 
about 86% would suffer tidal instability and merge into the single fast-rotating stars. 
This suggests that the dynamical stability limit for the observed W UMa systems is 
higher than the theoretical value, implying that the observed systems have probably 
suffered the loss of angular momentum due to gravitational wave radiation (GR) or 
magnetic stellar wind (MSW). 
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1 INTRODUCTION 



Interest in W UMa binaries was revived recently with the 
discovery of large number of new ones among blue strag- 
glers in open and globular clusters dKaluznv fe Sharall 19881; 
Kaluznvlll990l ; iMateo et al"1ll990l : lYan fc Matedll994lMateol 



1996I - ). It appears likely that at least some of blue stragglers 



are formed by the m erging of contact bin aries. The models 
of W UMa systems (|Li et alj |2004 120051 ) had shown that 
W UMa systems evolve into contact binaries with extreme 
mass ratios and then evolve into single, fast-rotating stars 
(FK Com stars) or blue stragglers due to Darwin's instabil- 
ity when the spin angular momentum of the syst em is more 
than a third of its orbital angular momentum (|Hutl [l980l : 
lEggleton fc Kiseleva-Eggletonll200lh . 

The spin angular momentum has a significant influence 
on the evolution of W UMa contact binaries, especially for 
those with extreme mass ratios. N eglecting the spin angular 
momentum of the small secondary, RasioJ (|l995T ) first derived 
a cutoff mass ratio for contact binary systems and the cutoff 
mass ratio depends on the structure and the contact degree. 
Meanwhile, iRasio fc Shapriol (|l995h identified the dynami- 
cal stability li mit at a conta ct degree of about 70% for W 
UMa systems. lYand l|2005bl ) derived a cutoff mass ratio of 
about 0.021 from the observed W UMa systems with mass 



ratios q < 0.25. In this work, using the observational data 
and the theoretical results reported recently, we investigate 
the dynamical stability limit of W UMa systems, and find 
that W UMa systems begin merging when their mass ratios 
are in a region of about 0.076-0.078. This suggests that W 
UMa systems with mass ratio q < 0.076 can not be observed. 
Meanwhile, we estimate that the observed W UMa systems 
with a minimum mass ratio of about 0.077, corresponding to 
a contact degree of about 86%, would suffer Darwin's insta- 
bility and begin merging. This suggests that the minimum 
mass ratio for the observed W UMa systems is larger than 
the theoretical value, implying that the observed W UMa 
systems have suffered the loss of angular momentum loss 
owing to GR or MSW. 



2 DYNAMICAL STABILITY LIMIT 

It is well known that a contact binary would suffer Darwin's 
instability and then coalesce into a fast-rotating single star 
(including FK Corn-type stars and blue stragglers) when the 
spin angular momentum of the system is more than a third of 
its orbital angular momentum. The spin angular momentum 
of a binary star can be expressed as 
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where Mi,2 and Ri t 2 are the masses and radii of the primary 
and the secondary in solar units, fci,2 the ratios of the gy- 
ration radii to the stellar radii for both components, and co s 
the spin angular velocity. The orbital angular momentum of 
the binary system reads 



J 01 



M 1 M 2 .2 

-A oj , 



Mi + M 2 



(2) 



where lj is the orbital angular velocity and A the orbital 
radius of the binary. We assume that the W UMa systems 
are in synchronous rotation (i.e. lu s = w„) and the ratios 
of gyration radii to the stellar radii for both components 
are equal (i.e. k\ — k\ = fc 2 ), the ratio of the spin angular 
momentum to the orbital angular momentum can be written 
as 



R 



J spin k 2 (MiR\ + M 2 R\ 



Jorb A 2 M 1 M 2 /(M 1 + M 2 ) 
, 2 l+q(Ri\ 2 f. . (R2^ 2 




(3) 



where q (= M2/M1 < 1) is the mass ratio of the binary 
system. If a contact binary system is a marginal contact 
one in which tw o components ha s just filled the inner Roche 
lobes, following lEggletonl (|l983T l. the relative radii of both 
components of the contact binary read 
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If the two components of a contact binary system have 
filled the outer Roche lobes, the relative radii of the 
components can be exp ressed by the following equations 
l|Yakut fe Eggletonll2005l ) 
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(6) 
(7) 



Using Eqs. (3), (4) and (5) and taking k 2 to be the values 
of main sequence stars (from 0.05 for no convection stars to 
0.21 for fully convective stars), we can obtain a relationship 
between the angular momentum ratio Ri and the mass ratio 
q for marginal contact systems 



Ri = k 2 
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setting equation Ri — 1/3, and solving (numerically) for 
q, we obtain the minimum mass ratio g m i n ,i n for stability 
of marginal contact binaries with the different values of k 2 
of MS stars and the relation between the minimum mass 
ratio g m in,in and k 2 is shown in Figure 1. Using Eqs. (3), 
(6) and (7), we obtain a relationship between the angular 
momentum ratio R and the mass ratio q for contact systems 
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Figure 1. The minimum mass ratio for stability g m i n in the 
Roche approximation. The parameter k 2 is the dimensionless gy- 
ration radius of the stars 



in which both components have filled the outer Roche lobes 
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setting equation R a = 1/3, and solving (numerically) for g, 
we obtain the minimum mass ratio, g m in,out, for stability of 
W UMa systems in which two components have filled the 
outer Roche lobes with the different values of k 2 and the re- 
lation between g m m,out and k 2 is also shown in Figure 1. As 
seen from Figure 1, the value of g m in,m is smaller than that 
of g m in,out at any value of fc 2 . This suggests that the stabil- 
ity of a contact binary with given masses indeed depends on 
the degree of contact /, or filling factor, < F = 1 — f < 1. 
Meanwhile, the more is the value of fc 2 , the larger are the 
minimum mass ratios (g m i n ,i n and g m i n ,out) and the differ- 
ence between g m i n ,in and gm in.out- 

We take fc 2 = 0.06 as iRasiol l| 19951 ). the ratio of the 
spin angular momentum to the orbital angular momentum 
against the mass ratio is plotted in Figure 2. Meanwhile, we 
obtain g m m,i n = 0.071 and g m i n ,out = 0.078 for fc 2 = 0.06. 
In general, the dynamical instability of W UMa contact bi- 
nary systems should be occurred at an overcontact stage 
with a higher degree of contact rather than at a marginal 
contact one. The numerical study of the equilibrium and 
stability of close binary systems l|Rasio fc Shapriol ll995T l 
had identified the dynamical stability limit at a contact 
degree of about 70%, corresponding to a minimum mass 
ratio of about 0.076. This suggests that the W UMa sys- 
tems with mass ratios in a region of 0.076-0.078 would 
suffer Darwin's instability and then coalesce into the fast- 
rotating single stars (including FK Com stars and blue strag- 
glers), implying that the W UMa systems with mass ra- 
tios q < 0.076 can not be observed. We collect the abso- 
lute parameters of some of W UMa systems (listed in Ta- 
ble 1) reported recently, tog et her with those com piled by 
iMaceroni fc van't Veer] ill996T ). lGazeas et ail l|2005h . and by 
lYakut & Egglctoi (|2005h 7 to determine the minimum mass 
ratio of the observed W UMa systems. Taking fc 2 = 0.06 
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Table 1. The physical parameters of some W UMa systems. 
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Figure 2. The angular momentum ratio R against the mass ratio, 
q. The solid lines represent theoretical results corresponding to 
the systems filling the inner Roche lobe and the outer Roche 
lobe, respectively. The solid dots represent the observed W UMa 
systems. 
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Figure 3. The angular momentum ratio R against the mass ra- 
tio, q. The solid line represents the fitting result. The solid dots 
represent the observed W UMa systems. 

as mentioned above, the observed systems are also plotted 
in Figure 2. It is seen in Figure that, except for two sys- 
tems with a larger scatter, all systems load in the region 
limited by the theoretical curves corresponding to the inner 
and outer Roche lobes, and none of them has surpassed the 
dynamical stability limit. However, it is difficult to obtain a 
direct polynomial relation between the ratio R and the mass 
ratio q by the least-squares fitting. But we have found a re- 
lation between the logarithms of the ratio R and the mass 
ratio q as the following 

logi? = -1.4347 + 0.099541ogg + 1.8506(log<j) 2 

+1.25287(log<j) 3 + 0.3291(logg) 4 , (10) 

the observed systems, together with the theoretical curve, 
are shown in Figure 3. As seen from Figure 3, the theo- 
retical curve has given a good fit to the observed systems. 
Setting equation logi? = log(l/3), and solving (numerically) 
for q, we obtain a minimum mass ratio, g m i n ,obs = 0.077, for 
the observed systems. It is larger than 0.076 estimated from 
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Refs. 




(days) 
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V410 Aur 


0.3663 


1.30 


0.146 


1.40 


0.61 


(1) 


XY Boo 


0.3705 


0.91 


0.186 


1.23 


0.61 


(1) 


AH Cnc 


0.3605 


1.21 


0.149 


1.36 


0.62 


(2) 


CW Cas 


0.3188 


1.06 


0.547 


1.01 


0.76 


(3) 


V776 Cas 


0.4404 


1.63 


0.13 


1.71 


0.71 


(1) 


V899 Her 


0.4212 


2.10 


0.566 


1.57 


1.22 


(5) 


BB Peg 


0.3615 


1.38 


0.362 


1.26 


0.76 


(6) 


V351 Peg 


0.5933 


1.63 


0.361 


1.87 


1.19 


(7) 


AU Ser 


0.3865 


0.90 


0.710 


1.10 


0.94 


(8) 


VZ Psc 


0.2613 


0.81 


0.800 


0.78 


0.70 


(9) 
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et al.(2005); (8) Gurol(2005); (9) Hrivnak et al.(1995) 

the theoretical result (|Rasio fc Shapriolll995l ). This suggests 
that the observed W UMa systems have probably suffered 
the loss of angular momentum. Since the loss of angular mo- 
mentum of the observed systems would result in the shrink- 
ing of their orbits and thus would cause their inner and outer 
Roche lobes to decrease and contact degree / to increase, so 
that leads / to have exceeded 70% before the observed sys- 
tems reaching the minimum mass ratio q = 0.076. 



3 DISCUSSIONS AND CONCLUSIONS 

The model s of low-mass W UMa-type systems (|Li et al.l 
|2004 I2005T ) have shown that W UMa systems would suf- 
fer Darwin's instability and evolve into fast-r otating single 
stars (i.e. FK Com stars or blue stragglers). iRasiol (|1995r ) 
predicted a cutoff of the mass ratio for W UMa systems 
at about 0.09. Considering the rotation of the secondary, 
we obtain that the minimum mass ratio for stability of the 
marginal contact binaries is of about g m i n in = 0.071, suggest- 
ing that the rotation of the secondary still has a definitely 
influence on the determi nation of the dynamical stability 
limit of W UMa systems. IRasiol (| 19951 ) had predicted that 
the dynamical stability limit should depend on the fill-factor 
of W UMa systems. Using the new results in studies of the W 
UMa systems, we investigate the dynamical stability limit of 
W UMa systems in which both components fill their outer 
Roche lobes and obtain a minimum mass ratio Qmin.out of 
about 0.078 for this kind of systems. This suggests that the 
dynamical stability limit of W UMa systems indeed depends 
on their fill-factor. 

The difference between g m i n ,m and gmin.out suggests 
that the dynamical stability of the W UMa contact bi- 
naries with given mass depends on the degree of contact 
[/ = (fi — Hm)/ (flo ut — fl m )1, wh i ch is consistent w i th the 
result predicted by IRasiol (|l995T ). iRasio fc Shapriol (|l995T l 
had identified the dynamical stability limit of W UMa sys- 
tems at about / = 70%, which corresponds to a mini- 
mum mass ratio of about 0.076. We obtain a minimum 
mass ratio of about 0.077 , corresponding to a contact de- 
gree / = (g min ,obs - gmin,m)/(<7mm,out - <7min,m) « 86%, for 

the observed systems. The difference between the Rasio & 
Shaprio's result and ours derived from the observed systems 
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suggests that the observed W UMa systems have proba- 
bly suffered the loss of angular momentum owing to GR or 
MSW during the evolution of the observed W UMa systems. 
In fact, Einstein soft X-ray observations and Inter national 
Ultraviolet Explorer (IUE) Ultraviolet observations ijVaianal 
Il98d : lEatonlll983l : iRucinski fe Vilhulll983l : IVilhulll983h have 
shown that W UMa systems are strong sources. This sug- 
gests surface activity of the kind observed on our Sun, and 
so presence of the magnetic field. Since the observed W UMa 
systems rapidly lose the angular momentum via MSW, the 
contact degrees of them have exceeded 70% before their mass 
ratios decrease to the minimum mass ratio q=0.076. 

Among the about 200 contact binaries with reliable pho- 
tometric data, the syste m V857 Her has an extremely small 
mass ratio, q = 0.0653 (Qian et al.l 120051 ) for observed W 
UMa systems at present. This is not only a new "record" 
for contact binary, surpassin g the well-known AW UMa with 
0.075 jRucinsklll992l ) and SV Crv with q = 0.066 



( Rucinski et al.l l200lh . It is a tremendous contribution for 
theory that W UMa systems with the extremely small mass 
ratios are found. Only a slight difference between observa- 
tions and theory not only indicates that theory have come 
close to a decent prediction, but also indicates that some 
physical processes taking place in the observed W UMa sys- 
tems have not be considered in the present theory. At first, 
the dynamical stability limit of W UMa systems may also 
depend on the structure (i.e. the value of fc 2 ) and the value 
of fc 2 may decrease with the evolution. In addition, the dif- 
ferential rotation of the observed systems would lead the 
ratio of the gyration radii to the stellar radii of the compo- 
nents of the W UMa systems to be decreased. Therefore, the 
differential rotation of the components of the W UMa sys- 
tems is another physical reason why some W UMa systems 
with the mass ratios smaller than the minimum mass ratio 
0.076 still can be observed. However, contact binaries with 
mass ratios significantly below that of SX Crv and V857 Her 
should not be observed, since they are tidally unstable and 
quickly merge into a single, fa st-rotating o bject, on a tidal 
timescale of about 10 3 -10 4 yr (lRasiolll995l ). 
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